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Abstract 



Using data collected in the region of the T(4S) resonance with the CLEO 
II detector operating at the Cornell Electron Storage Ring (CESR), we present 
the first observation of B mesons decaying into the charmed strange baryons 
~° and E+. We find 79 ± 27 E° c and 125 ± 28 H+ candidates from B decays, 
leading to product branching fractions of B(B — > 'B^X)B(E^ — > H~7r + ) = 
(0.144±0.048±0.021) x 10" 3 and B(B -» Z+X)B(Z+ -» H-7r+vr+) = (0.453± 
0.096 i°;°|)xl0- 3 . 
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Charmed baryon production from the decays of B mesons has been previously reported 
by ARGUS QXJ] and CLEO |2]|| . Here, we report the first observation of the charmed-strange 
baryons S° and from B decays M, which have previously been observed only in direct 
charm production |5HT0 . 



In e + e~ annihilations at the T(4S) resonance (10.58 GeV), charmed baryons can be 
produced either from B meson decay or from hadronization of cc quarks produced in the 
continuum. Since the b quark couples predominantly to the c quark, B meson decays to 
the charmed strange baryons S° (csd) and 5+ (csu) will proceed through either spectator or 
exchange diagrams. Decays mediated by the coupling b — ► cW~ with W~ — > ud produce final 
states of the form E c YX h and E C NX S , where Y is a hyperon (A, E, H, etc.), N is a nucleon, 
and X h (X s ) denotes non-strange (strange) multi-body mesonic states (see Figure 0(a)). As 
shown in Figure Mb), decays mediated by b — > cW~ with W~ — > cs can lead to states of 
the form S c @ c ||iT| , [l2| , where G c denotes any charmed non-strange baryon. The authors 
of Refs. [T3| and [IT4I] predict branching ratios of (1.0 — 1.8) x 10~ 3 for those decays. The 



process b — > uW~ with W~ — > cs leads to final states of the form S C F, but should be highly 
suppressed by the small b — > u coupling. 

There are several theoretical calculations that attempt to derive the two-body contri- 
bution to charmed baryon production in B decays. In the diquark model fl3| baryons of 
spin |(|) are modeled as bound states of quarks and scalar (vector) diquarks. The b quark 
decays to a scalar diquark and an antiquark; the latter combines with the light antiquark 
accompanying the b quark to form an antidiquark. The creation of a qq pair then leads to a 
baryon and antibaryon in the final state. The authors of Ref. |14| calculate decay amplitudes 
based on QCD sum rules, replacing both the B meson and the charmed baryon in the final 
state by suitable interpolating currents. There are also treatments that determine the rates 



for exclusive baryonic B decays in terms of three reduced matrix elements |L5j , on the basis 



of the quark diagram scheme [flB] , using the constituent quark model |]17|1 , and using the 
pole model ||18[ . The latter four calculations do not quote explicit predictions for branching 
fractions of B decay modes which yield S c baryons. 

For this analysis we used 3.1 fb _1 of data taken on the T(4S) resonance, corresponding 
to 3.3 million BB events. To estimate and subtract continuum background, 1.6 fb _1 of 
data were collected 60 MeV below the resonance. The data were collected with the CLEO 
II detector operating at the Cornell Electron Storage Ring, CESR. The CLEO II detector 
|19| is a general purpose solenoidal-magnet detector with excellent charged particle and 
shower energy detection capabilities. The detector consists of a charged particle tracking 
system surrounded by a scintillation counter time-of-flight system and an electromagnetic 
shower detector consisting of 7800 thallium-doped cesium iodide crystals. These detectors 
are installed within a 1.5 T superconducting solenoidal magnet. Incorporated in the return 
yoke of the magnet are chambers for muon detection. 

Charge measurements from the drift chamber wires provide specific ionization loss 
(dE/dx) information. To obtain hadron identification, dE/dx and available time-of- 
flight (TOF) measurements are combined to define a joint x\ = [{(dE/dx) meas — 
(dE/dx) cxp }/a dE/dx ] 2 i + [{(T) meas - (T) cxp } /a TOF ]f , where i corresponds to the pion, kaon, 
and proton hypotheses. A x 2 -probability is then calculated for each hypothesis, and par- 
ticle identification levels for each of the hypotheses are derived by normalizing to the sum 
of the three probabilities. A particle is identified with a specific hypothesis if its particle 
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identification level for it is greater than 0.05. 

We reconstruct S° (H+) candidates through the decay chain S° — > (5+ — > 

S _ 7r + 7r + ), H~ — > Att~ , and A — > £>7r~. We study the S c momentum spectra using the 
scaled momentum x p = p/(E^ eam — m| ) 1//2 , where p and m= c are the S c momentum and 
mass, respectively, and -Ebeam is the beam energy. We require x p < 0.5, the kinematic limit 
for H c baryons produced from B decays. This requirement reduces the background from 
continuum cc. 

The A candidates are formed from pairs of oppositely charged tracks, assuming the higher 
momentum track to be a proton and the lower momentum track to be a pion. We also require 
the higher momentum track to be consistent with the proton hypothesis. The invariant mass 
of A candidates has to be within 5.0 MeV/c 2 (corresponding to 2.5 standard deviations) of 
the known A mass. We have not required A candidates to point towards the primary vertex, 
since A's decaying from H~'s can travel as much as a few centimeters before decaying and 
can have appreciable impact parameters. To reduce the background from tracks coming from 
the interaction point, we require the radial distance of the A decay vertex from the beam 
line to be greater than 2 mm. 

The S~ candidates are formed by combining each A candidate with the remaining neg- 
atively charged tracks in the event, assuming the additional track to be a pion. The decay 
vertex of the S~ candidate is reconstructed by intersecting the extrapolated A path with the 
negatively charged track. We require the radial distance of the S~ decay vertex from the 
beam line to be greater than 2 mm and less than the radial distance of the A decay vertex. 
In addition, the reconstructed S~ momentum vector has to point back to the interaction 
point. The invariant mass of the S~ candidates has to be within 6.5 MeV/c 2 (corresponding 
to 3 standard deviations) of the known S~ mass. 

To reconstruct S° candidates, we form combinations of H~ with one positively charged 
track, and to reconstruct candidates, we combine each H~ with two positively charged 
tracks. These additional charged tracks are required to originate from the interaction point 
and to be consistent with the pion hypothesis. 

To find the S c signal yields, we fit each invariant mass distribution to the sum of a Gaus- 
sian function of fixed width and a second order polynomial background, both for the T(45) 
and the continuum data. The fixed widths for the two modes were determined using a Monte 
Carlo simulation of the detector, resulting in widths of 8.0 and 6.8 MeV for the H° and the 
S+, respectively. We scale the continuum yields to account for the differences in luminosities 
and cross sections in the two data sets with the scale factor (£Y(4S)/£cont)(^ra n t/-^T(4S))> 
where £~r(4S) an d C cont are the luminosities, and £7-r(4S) an d E cont are the beam energies on 
the T(45*) and on the continuum. Figure shows the invariant mass distributions of the 
2~7r + and H~7r + 7r + combinations from T(45) and scaled continuum data. After subtracting 
the scaled continuum yield from the T(AS) yield, we observe 79 ± 27 S° candidates and 
125 ± 28 5+ candidates from B decays. The errors are statistical only. The fitted S c masses 
are consistent with the current world averages. 

To measure the product branching fractions for the two decay modes, we divide both data 
and Monte Carlo into x p intervals. The reconstruction efficiency in each mode is found as a 
function of x p using Monte Carlo simulations. Tables | and [IJ show the continuum subtracted 
raw yields y r (x p ) and efficiency-corrected yields y c (x p ). We also give the fractional decay rate 
in each x p interval, (l/NB)( y dy c /dx p ), where Nb is 2N B - B - ) for S° and production. We find 
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B(B -> E° C X)B(E° C -> S-7T+) = (0.144 ± 0.048 ± 0.021) x 10~ 3 and -6(5 -> 



7T 



7r + ) = (0.453 ±0.096 1 oioes) x 10 , with the first error being statistical and the second 
being systematic. The main sources of systematic error are due to uncertainties in the 
reconstruction efficiencies for A (5%) and S~ (7%), variations in the selection criteria (8-9%), 
uncertainties in particle identification (5%), charged particle tracking (1% per track), and 
the Monte Carlo predictions for the signal width (4%). These result in a total systematic 
uncertainty of about 14%. In addition, we assign a +12% systematic uncertainty in the 
2~7r + 7r + case for the possible resonant substructure S*°7r + , since this would decrease the S+ 
reconstruction efficiency considerably. 

We can convert these product branching fractions into absolute branching ratios using 
the following branching fractions of S° — > 2~7r + and — > S _ 7r + 7r + , derived by CLEO 
0: B(~° c -> H-7T+) = / 5L / Hc (0.52 ± 0.16 + _° f )% and B(E+ -> = / 5 l/h c (2.5 ± 

0-6 loJ)^) where fs c = B(E C — > H£ + z/j)/£>(H c — > £ + X) < 1 (current predictions range from 
0.4 to 0.9 0,|22|), and / Si = (r|| / /r^)(r^/rf L ), with T 5L being the total semileptonic 
width. These numbers are actually slightly different from the published values, since we are 



now using an updated value for rf L = (0.165 ± 0.009) ps~ x f23|j24[ (instead of the previous 
value of (0.138 ± 0.006) ps _1 ). In addition, we have introduced the factor fsL to account 
for predictions of the semileptonic width of the S c being quite different from that of the 
A c |25j (2 to 3 times as large), which in turn should be different from that of the D [26 



namely about 1.5 times as large. This leads to the following absolute branching ratios: 
B(B -> ~° C X) = f s lf Ec \2.8 ± 1.3+8;?)% and B(B - H c +X) = f s lf~ c \l.8 ± 0.6 t °: 4 7 )%. 

In Figure ||| we present the corresponding efficiency-corrected momentum spectra of S° 
and 5+ baryons in B decays. Superimposed on the measured spectra are the results from 
Monte Carlo simulations of the decays B — > H c A( c )(n7r), n — 0, ...,3. Comparing the mea- 
sured spectra with Monte Carlo predictions indicates that two-body final states such as S C A 
and S C E are suppressed while multi-body final states seem to be dominant. We are not yet 
sensitive to b —>■ ccs decays leading to final states of the form S C A C or S C S C , which are pre- 
dicted by the authors of Refs. |13| and |14] to have branching fractions of only (1.0—1.8) x 10~ 3 
for those decays. These branching fractions are at least an order of magnitude lower than 
the inclusive branching fractions for B — > S C X. 

In summary, we have presented the first observation of B mesons decaying into the 
charmed strange baryons S° and 5 + . From an examination of the measured H° and 
momentum spectra, it is not clear which of the possible production mechanisms b — > cud 
or b — > ccs is preferred or dominant, since the observed momentum spectra are consistent 
with both mechanisms. It seems, however, that decays involving a heavier anti-baryon or 
multi-body decays are favored. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftung, Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, and the A. P. Sloan Foundation. 
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FIG. 1. Possible B — ► baryon decay mechanisms: (a) B — > c iVX and S c yi , (b) i? — > H c c X 
and B — > yH c X; iV stands for any non-strange non-charmed baryon, V for any strange and 
non-charmed baryon, and C for any charmed and non-strange baryon. 
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FIG. 3. Efficiency-corrected momentum spectra for (a) S° and (b) from B decays. The 
superimposed curves indicate the spectra derived from Monte Carlo simulation of the decays 
B — > H c A( c )(n7r), n = 0, ...,3. The Monte Carlo curves have been normalized to data, except 
for the two-body decays, where the normalization is arbitrary. 
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TABLE I. Inclusive H° production in B decays. 



Ax p 


Raw yield 


Corr. yield 


(l/N B )(dy c /dx p ) 




Vr(Xp) 


Vc{x p ) 


[io- 3 ] 


0.0 - 0.1 


27.0 ±6.5 


358.8 ±88.1 


0.54 ±0.13 


0.1-0.2 


33.4 ± 13.5 


399.5 ± 162.3 


0.60 ± 0.24 


0.2-0.3 


43.5 ± 13.6 


482.8 ± 152.5 


0.72 ±0.23 


0.3 - 0.4 


-18.1 ± 12.2 


-191.5 ± 129.5 


-0.29 ±0.19 


0.4 - 0.5 


-6.9 ±13.3 


-89.7 ±174.1 


-0.13 ±0.26 


0.0 - 0.5 


78.9 ± 27.2 


959.9 ±323.1 





TABLE II. Inclusive £+ production in B decays. 



Ax p 


Raw yield 


Corr. yield 


(l/N B )(dy c /dx p ) 






Vc(Xp) 


[IO" 3 ] 


0.0-0.1 


10.0 ±7.0 


417.1 ± 295.0 


0.62 ± 0.44 


0.1 - 0.2 


47.0 ± 14.3 


1273.5 ± 392.6 


1.91 ±0.59 


0.2 - 0.3 


41.8 ±13.0 


901.4 ± 285.5 


1.35 ±0.43 


0.3-0.4 


20.2 ± 13.6 


344.2 ± 232.8 


0.52 ±0.35 


0.4-0.5 


6.0 ± 12.4 


89.6 ± 186.0 


0.13 ±0.28 


0.0-0.5 


125.0 ±27.6 


3025.8 ± 641.5 
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